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ABSTRACT 

We show that r-mode oscillations distort the magnetic fields of neutron stars and that their occurrence is likely 
to be limited by this interaction. If the field is > 10^^ (il/ils) G, where fl and Hb are the angular velocities of 
the star and at which mass shedding occurs, r-mode oscillations cannot occur Much weaker fields will prevent 
gravitational radiation from exciting r-mode oscillations or damp them on a relatively short timescale by extracting 
energy from the modes faster than gravitational wave emission can pump energy into them. For example, a 10^" G 
poloidal magnetic field that threads the star's superconducting core is likely to prevent the £ — 2 mode from being 
excited unless fl exceeds 0.35 fis. If fl is larger than 0.35 ils initially, the £ = 2 mode may be excited but is 
likely to decay rapidly once fl falls below 0.35 fis, which happens in ^ 15"^ if the saturation amplitude is ^ 0.1. 
The r-mode oscillations may play an important role in determining the structure of neutron star magnetic fields. 

Subject headings: stars: neutron — stars: oscillations — instabilities — stars: rotation — magnetic fields 



1. INTRODUCTION 

The properties of r-mode oscillations in rotating Newto- 
nian stars were initially investigated by Papaloizou & Pringle 
(1978), who suggested they may explain the short-period os- 
cillations seen in cataclysmic variables. Recent work demon- 
strating that gravitational radiation can cause the r modes of 
rotating fluid stars to grow (see, e.g., Andersson 1998; Anders- 
son, Kokkotas, & Stergioulas 1998; Friedman & Morsink 1998; 
Lindblom, Owen, & Morsink 1998) has generated renewed in- 
terest in these modes (see Bildsten 1998; Madsen 1998; Owen 
et al. 1998; Andersson, Kokkotas, & Schutz 1999; Kokkotas 
& Stergioulas 1999; Levin 1999; Lindblom & Mendell 1999; 
Lindblom, Mendell, & Owen 1999; Kojima 1998; Kojima & 
Hosonuma 1999). Unlike nonaxisymmetric density perturba- 
tions, the r modes would be unstable at an arbitrarily small 
angular velocity if neutron stars were made of perfect fluids. 
This increases the potential astrophysical significance of the 
r modes. The studies carried out to date have established the 
basic features of the r-mode instability in fluid stars, but im- 
portant questions remain. These include interaction of r mode 
oscillations with the stellar magnetic field and the proton flux 
tubes in the core; the effects of core-crust coupling; and the 
nonlinear development of these oscillations. 

Here we investigate the distortion of the stellar magnetic field 
caused by r mode oscillations and the effect of this distortion 
on the onset and evolution of the r-mode instability when the 
effects of the superconductivity of the neutron star core are in- 
cluded. We show that a magnetic field of ~ 10^^ G will prevent 
r-mode oscillations from occurring, even if the star is spinning 
very rapidly. A much weaker field will prevent r-mode oscil- 
lations from being excited by gravitational wave emission or 
cause them to die out as the star spins down. A more detailed 
discussion of the fluid motions produced by the r modes and 
their interaction with the magnetic field of the star will be pre- 
sented in a longer paper (Rezzolla et al. 2000). 
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2. PROPERTIES OF H-MODE OSCILLATIONS 

The lowest-order Newtonian velocity perturbations produced 
by small-amplitude r-mode oscillations of a slowly and uni- 
formly rotating, nonmagnetic, inviscid, fluid star can be deter- 
mined by expanding the fluid equations in powers of the di- 
mensionless mode amplitude a and the ratio ^I/Ub, keeping 
only the lowest-order nonvanishing terms, which are of order a 
and (il/ris)^. Here is the angular velocity of the star and 
Qb = (2/3)(7rG'p)^/^ is approximately the angular velocity at 
which mass-shedding is expected to occur (see Lindblom et al. 
1998). The Eulerian velocity perturbations that are solutions of 
these linearized fluid equations for barytropic stars are of axial 
type and may be written 6vi (r, 6, 4>, t) = a^lR {r/Rf Y^e"'^, 

where R is the radius of the star, Y^^^ is the magnetic-type 
vector spherical harmonic, and a is the angular frequency of 
the mode in the inertial frame; to this same order, 5v'[ = 
(see Provost, Berthomieu, & Rocca 1 98 1; Friedman & Morsink 
1998; Lindblom etal. 1998). 

We are interested in the velocities of fluid elements in the 
frame corotating with the unperturbed fluid of the star. These 
velocities can be determined by using 5vi to compute the mo- 
tion of a given element of fluid. In a spherical coordinate basis 
(r, 9^ (p), the velocity of the fluid element located at r, 6, cf) is 

e{t,e,(j)) = aujf{r)ciSm0{sm0Y"^cos{£(l3 + ujt) , (1) 

<i){t, e, (f)) = -aujf{r)ce cos 6i(sin ey-"^ sin {£(j) + ujt) , (2) 

where the dot denotes the total derivative with respect to 
time, /(r) = {r/RY^^ is the radial mode function, q = 
(-l/2)^[(^ - iy]-^[{2£ + l)!/47r]i/2, and = a + ^ = 
2f2/(£ + 1) is the angular frequency in the corotating frame. 

We expect r waves to cause fluid elements to drift as well as 
gyrate, just as sound and water waves cause fluid elements to 
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drift as well as oscillate (see Landau & Lifshitz 1959, § 64; the 
time average of 6vi at sl fixed location in the corotating frame is 
of course zero). We expect the drift of fluid elements given by 
the velocity field 6vi to be qualitatively correct and perhaps ex- 
act to order a^. (Computing the drift of fluid elements caused 
by sound waves and shallow water waves using the solutions 
Svi of the relevant linearized fluid equations as we have for the 
r waves gives these drifts exactly to second order in the wave 
amplitude; see Rezzolla et al. 2000.) As Figure 1 shows, the 
and 9 excursions of a fluid element during a gyration are compa- 
rable. The drift velocity Vd produced by 6vi is in the azimuthal 
direction; its sign and magnitude vary with r and 9. Fluid el- 
ements in the northern and southern hemispheres equidistant 
from the rotation equator drift in the same direction at the same 
rate, although their gyrational motions are different. 

An analytical expression for the drift velocity of a given fluid 
element can be obtained by expanding equations (jl|) and (^ 
in powers of a, averaging over a gyration, and retaining only 
the lowest-order nonvanishing term. For the £ — 2 mode, this 
gives Vd{r,9,t) = H2{9)a'^{t)uj{t)R{r / Rf4>, where K2{9) = 
(l/2)''(5!/7r)(sin^ 9 - 2 008^ 9) and ^ is the unit vector in the 
direction. Hence the net displacement in the azimuthal direc- 
tion after an oscillation is approximately 27ra times the radius 
of gyration. The total displacement in from the onset of the 
oscillation at to to time t is 

A0(r,0,i) = k{9) (^^yj\\t')uj{t')dt' . (3) 

Direct numerical integrations of equations (|l]) and (^ show that 
expression (|]) accurately describes the secular displacement 
given by these equations, even for a = 1. 

3. DISTORTION OF THE STELLAR MAGENTIC FIELD 

The r modes produce fluid motions perpendicular to any 
plausible stellar magnetic field. Such motions distort the field, 
increasing its energy. As we show in § 4, the back-reaction of 
the magnetic field on the r modes can be related to the increase 
in magnetic field energy. At the core temperatures <^ 10^ K 
relevant to r-mode excitation in newly formed neutron stars, 
the electrical conductivity of the core is ^ 10^^ s^^ (see Lamb 
1991, § 4.2), so the diffusion time for fields that vary on length 
scales > 10^ cm is ^ 10^ yr, much longer than the timescales of 
interest here. The diffusion time for older neutron stars is much 
longer. We therefore neglect diffusion and compute the effect 
of the fluid motions on the magnetic field in the ideal MHD 
limit in which field lines are "frozen into" the fluid and carried 
with it. The evolution of B can be computed to leading order 
in a without considering the change in p, because V • Svi — 0. 
Here we also neglect any instability of the magnetic field. We 
defer consideration of the back-reaction of the magnetic field 
on the oscillations to § 4. With these approximations, the evo- 
lution of the magnetic field is given by the induction equation 
dB/dt = W X {ux B),w\\er& d/ dt = {d/dt + dxr-V ~nx) 
describes the rate of change of S at a fixed position in the frame 
corotating with the star and u is the fluid velocity in this frame. 

The distortion of the magnetic field during an 7--mode gy- 
ration is given by the induction equation with u = 6vi. The 
perturbation of the field during a single oscillation is 6B « 
{dB/dt){TT/n), where [dB/dt) w 5viB/R w a^B, because 
Svi has magnitude ailR and varies with radius and latitude on 
the scale R. Hence the distortion of the field during an oscilla- 
tion is SB w naB. 
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Fig. 1 . — Motions of eight fiducial fluid elements in the northern hemisphere 
of a star undergoing £ = m = 2 oscillations with amphtude a = 0.1. The 
projected trajectories sin6'(t) cos and sin 9(t) cos <?!>(t) show 
how the motions over five oscillation periods 2tt/u! would appear to a distant, 
corotating observer in the rotation equator of the star. The initial positions 
of the fluid elements are indicated by the 'star' symbols; their positions af- 
ter each oscillation period ai'e indicaled by the filled squares. The trajectories 
were computed by integrating eqs. (jlh and (fcl) numerically. 

The magnetic field of the star is increasingly distorted by r- 
mode oscillations, because Vd varies with radius and latitude. 
The shear produced by Vd generates an azimuthal magnetic field 
from any poloidal field Bq that is initially present. The evolu- 
tion of the azimuthal field is given by the induction equation 
with u ~ Vd ~ Ka^uRcj), so the change in the azimuthal field 
from the time to that the oscillation begins to time t is given ap- 
proximately by AB"^ (t) « J^^ Ka^{t)Lo{t)B^dt. The change in 

B'^ can be computed more accurately by noting that the induc- 
tion equation can also be written d'B/dt = (B ■ W)Sv, where 
d' /dt = d/dt + 6v ■ W is the rate of change of B in the frame 
comoving with the fluid. This equation can be integrated to re- 
late the magnetic field at the position a; of a given fluid element 
at time t to the magnetic field at the position xo of the same 
fluid element at the earUer time to (see Parker 1979, § 4.3). The 
result is 

B=(x,t) = B\xoM)^^, (4) 

where dx^ /Oxq is the coordinate strain between to and t. 

The evolution of the azimuthal magnetic field can be com- 
puted by averaging the azimuthal component of equation (Q) 
over a time At much longer than a single gyration but much 
shorter than the time for a fluid element to drift around the star. 
Symmetry arguments and direct numerical integration of equa- 
tions (|l]) and (^ show that the gyration-averaged displacement 
is independent of (j!) and hence that 9(a;'^(i))/9(}!)(io) = 1- Thus 
the gyration-averaged azimuthal field at position x at time t is 

{B^ (f , t)) = Bi + Bl ^^AiZ + Bl , (5) 
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where Bq = B^XQ,to) and the brackets indicate that the 
quantity is gyration-averaged. Hereafter we will consider only 
gyration-averaged quantities, omitting the brackets for brevity. 

As explained in § 4, the damping effect of the magnetic field 
can be related to the rate at which the energy in the magnetic 
field increases. After a short time, the energy in the magnetic 
field is dominated by the volume average of \B'^{x, which 
depends on the structure of the initial magnetic field as well as 
on the r mode in question. \B'^{x, t) \ can be computed numer- 
ically from equation (0) or analytically from equation using 
equation (^) with k{9) replaced by its appropriately weighted 
latitude average R. For the £ = 2 mode and an initial magnetic 
field that is dipolar and aligned with the rotation axis of the star, 
K — 0.3; for other modes and field orientations, R is different 
but of the same order Hereafter, in computing volume averages 
we assume a dipolar initial field and exclude the volume inside 
r = -^i? because of the divergence of the dipolar field as r ^ 0. 

We illustrate the possible effects of the interaction of r modes 
with the magnetic field of a newly formed neutron star by 
considering the evolutionary scenario outlined by Owen et 
al. (1998). In this scenario the star is assumed to be completely 
fluid, the initial spin rate is assumed to be Hb, and the sat- 
uration amplitudes of the r modes are assumed to be « 1. The 
velocity field is taken to be Svi, even though this expression 
is valid only for ft^ <^ and a <C 1. Owen et al. assumed 
that viscous damping is small enough initially that gravity wave 
emission can cause velocity perturbations of r-mode character 
to grow, but that gravity wave emission has no other effect on 
the velocity field. The growth timescale is assumed to be the 
gravity wave emission timescale ItgwI^ which is several tens 
of seconds for the fastest-growing, £ = 2 mode. Owen et al. 
assumed further that exponential growth ends when the ampli- 
tude of this mode reaches its saturation value, that this is the 
only nonlinear effect, and that the star then spins down at the 
rate given by angular momentum loss to £ = 2 gravitational 
radiation. In their scenario, viscous damping of the r modes 
becomes important after about a year 

We have calculated the evolution of the stellar magnetic field 
in this scenario, except that for simplicity we have neglected 
viscosity. Although we mention here our result for = ils, 
in order to allow comparison with the results of Owen et al. 
(1998), we emphasize that the approximations made in deriving 
Svi are accurate only for f2 -C f^s . A poly tropic stellar model 
wifli r = 2 and mass M = 1.4 A/q has R = 12.5 km and 
flB = 5.60xl0-'^rads-i = n*^. Forf^o = the ^ = 2 mode 
initially grows exponentially on the timescale ItqwI = 37 s. If 
the amplitude of the initial perturbation is 10~^ and the satu- 
ration amplitude is 0.1, the mode saturates after about 430 s. 

The volume-averaged azimuthal magnetic fields (AB'^) at sat- 
uration and after 1 year are - 300 (B^) G and - 10^ (B^) G, 
respectively, where the hats indicate an orthonormal basis and 
(Bq) is the volume-averaged poloidal component of the initial 

magnetic field. The evolution of (AS*) computed analytically 
using equations ^ and (^ agrees well with the evolution com- 
puted numerically using equations ([l|), (^, and (Q). 

4. EFFECT OF THE MAGNETIC FIELD ON R-MODE OSCILLATIONS 

The evolution of 7--mode oscillations in magnetic stars is 
likely to be affected not only by gravitational radiation and vis- 
cous dissipation, but also by the loss of rotational and mode 
energy to electromagnetic radiation and the growth of the stel- 
lar magnetic field. Here we focus on the back-reaction of the 



magnetic field on the £ = 2 mode. Although a detailed analy- 
sis of the change in the character and evolution of the r modes 
caused by their interaction with the stellar magnetic field is be- 
yond the scope of the present paper, the important qualitative 
effects can be estimated simply. If the magnetic field is strong 
enough, it will prevent r-mode oscillations from occurring. A 
much weaker magnetic field will prevent gravitational radia- 
tion from exciting r-mode oscillations or cause them to die out 
as the star spins down. The estimates presented here are con- 
firmed by more detailed calculations (Rezzolla et al. 2000). 

Prevention of r-mode oscillations. — The distortion of the 
magnetic field during an r-mode oscillation will prevent the os- 
cillation from occurring if the energy SEm that it costs to distort 
the magnetic field in the manner required by the mode is greater 
than the energy E in the mode. (The fractional change in the 
mode energy produced by gravitational wave emission during a 
single oscillation period is 4tt/uj\tgw\- It is always ^ 1 and 
therefore can be neglected in this comparison.) We denote the 
critical magnetic field that prevents r-mode oscillations i?crit,p- 

Consider first 5Em for a normal stellar core (we discuss 
the more realistic but less familiar case of a superconduct- 
ing core below). An estimate of the change in the magnetic 
energy of the star during an oscillation is Jy{6B^ /8tt) dV, 
where 6B is the magnetic field perturbation. As shown in § 3, 
SB w naB and hence SE'^^ « (iraB)^ /30. The mode en- 
ergy E is Jy{p6vi/8n) dV (see Lindblom et al. 1998), where 
Svi « am{r/RY, and hence E « {?,/28'n:){anRfM. Thus 
SE"^ >E If B> Bi^it.p « 4 X IQie {n/n*s)Ml'] R-^'^ G; 
hereMi.4 = M/1.4MQand i?i2.5 = i?/12.5km. 

The core of the neutron star is expected to be a Type II su- 
perconductor when the temperature there falls below the crit- 
ical temperature - 3 x 10^ K (see Lamb 1991, § 4.1). 
This is expected to occur after a few days (see Tsuruta & 
Nomoto 1988). In the superconducting core, the magnetic field 
is confined to flux tubes. The total energy of the flux tubes is 
{(j)totHf/8TT)£f, where 0tot ~ ttR^B^ is the total magnetic 
flux, Hf is the field inside a flux tube, and ^/ is the average 
length of a flux tube. During an oscillation the flux tubes are 
sheared a distance « dvi{TT/uj), which increases their length 
by 6£ « {l/2){37T/2)^a'^(r/R)^R and their total energy by 
SE^f' = {l/7){3Tra/8fBHfR^. If B < « 10^^ G, 
then Hf « Hcu but as S ^ Hc2 ~ 10^'' G, the flux fiibes 
crowd one another and Hf Hc2- When B = Hc2, the 
normal cores of the flux tubes touch and the superconducting 
state is destroyed (see de Gennes 1966, Chs. 3 & 6). If fl 
is « fls, the magnetic field needed to prevent oscillation of 
the £ ^ 2 mode is - lO^^ G and hence B^rit.p ^c^it.p^ 
whether or not the core is superconducting. If instead il <C tls, 
thenB^^gp « 10^^ {n/0.0in*B)^H-^\^Mi,4Ri2,5G, where 
i/,i,i6 = i^ci/10i6G. 

Damping of r-mode oscillations. — The increase in the en- 
ergy of the magnetic field produced by the r-mode drift re- 
duces the mode energy. If the rate dEm/dt at which the mag- 
netic field extracts energy from a mode is greater than the rate 
dEcR/dt at which gravity wave emission pumps energy into 
it, the mode either will not be excited or, if it was excited ear- 
lier, will decay. For a given stellar magnetic field, the condition 
dEoR/dt — dEm/dt defines a critical spin rate fJcrit above 
which r modes are excited and below which they are damped. 

If the core were normal, dE^/dt would be « 

{l/4Tr)JyB^{dB^/dt)dV, once B^ exceeds BP. As 
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shown in § 3, {dB^/dt) « Ka^Lj{r/R)^BP, so dE^/dt « 
{K/45)a'^il,{B't'BP)R^ and the critical spin rate would 
be l^^it « 0.31]Jj(B'>BP/1026G2)i/7M-4'/^ii'^2^/^ in- 
dependent of the mode amplitude. In reality, the core 
is expected to be superconducting and dEf^ /dt « 

{(j)totHf/8TT){def/dt) « {(j)t,tHf/8TT)va, once B^ ex- 
ceeds BP. Hence the actual critical spin rate is likely to be 

0.7 n*s {BP /W^G)^/^ {Hj /WGy/Hl-l'' R-l'^\ again 
independent of the mode amplitude. 

5. DISCUSSION 

The results reported here indicate that interaction of r-mode 
oscillations with the stellar magnetic field is likely to be impor- 
tant and may impose significant constraints on the occurrence 
and persistence of these oscillations in magnetic neutron stars. 
Our analysis shows that if the poloidal magnetic field of the 
star is ^ 10^*^ G, the increase in magnetic energy required to 
keep the 1 = 2 mode alive will exceed the energy that gravity 
wave emission can supply, unless the stellar angular velocity 
f2 exceeds 0.35 Jig. For the model of a newly formed neutron 
star discussed in § 3, falls below 0.35 ils in <^ 15*^ if the 
saturation amplitude is ^ O.L Once damping of an existing 
r mode begins, the mode decays quickly because gravity-wave 
emission decreases rapidly as the spin rate continues to decline. 

The energy density in a magnetic field of the strength needed 
to damp the r modes is 10^® of the energy density in the 
mode (the smallness of this number reflects the weakness of 
pumping by gravity wave emission). In contrast, the ratio of 
the fluid displacement after one gyration to the radius of gy- 
ration is {2TTVd/oj)/ {aR) « 2TrKa sa 6a, which is many or- 
ders of magnitude larger for mode amplitudes of interest. It 
therefore seems unlikely that the magnetic field will distort the 
r mode velocity field enough to cancel the fluid drift before 
it suppresses the oscillation. If r modes grow to appreciable 
amplitudes before the crust solidifies, their interaction with the 
initial poloidal magnetic field may generate strong azimuthal 
fields in young neutron stars. 

Spruit (1998) assumed that angular momentum loss to r- 
mode gravitational radiation would eventually cause the interi- 
ors of old, accreting neutron stars with millisecond spin periods 
to rotate half as fast as their outer layers. Neglecting the solid 
crust and its coupling to the interior, he estimated that this dif- 
ferential rotation would generate an azimuthal magnetic field 



^ 10^^ G before growth of the field ends. The results reported 
here indicate that if r mode oscillations occur, then the fluid 
motions they produce will themselves generate an azimuthal 
magnetic field. However, the damping effect of a moderately 
strong poloidal field will prevent such oscillations from devel- 
oping. 

Important issues remain unresolved. One is whether the 
back-reaction of the gravitational radiation force affects only 
the amplitudes of the oscillations or also the structure of the 
velocity field in the star We note that if it is rotating suffi- 
ciently rapidly, a uniformly rotating, incompressible (Maclau- 
rin) spheroid subject to a nonaxisymmetric perturbation is 
driven by (mass quadrupole) gravitational radiation to a non- 
radiating, triaxial configuration of uniform vorticity (a "station- 
ary football with internal mass currents"; see Miller 1974; De- 
tweiler & Lindblom 1977; Lai & Shapiro 1995). 

Another important question is the nonlinear development and 
saturation of the r modes. Do the low-^ modes grow to ampli- 
tudes of order 1 , or do they stop growing earlier because of the 
loss of mode energy to higher-^ modes? We note that pump- 
ing of r modes by gravitational radiation is a relatively weak 
effect, so that even a small energy loss rate due to mode-mode 
coupling can balance it. 

As the magnetic field grows, instabilities may allow mag- 
netic flux to leave the star, disrupting the r modes and limiting 
the strength of the azimuthal field. This is most likely in hot, 
young neutron stars in which relaxation to (3 equilibrium is rel- 
atively fast and the effects of stable stratification correspond- 
ingly reduced. Other important issues include the effects on the 
r modes of the solid crust, which is expected to begin forming 
when the surface temperature falls below ^ 3 x 10^° K (see 
Lamb 1991, § 3), and the interaction of the neutron vortices 
in the core with the proton flux tubes and the crust (see Sauls 
1989; Lamb 1991; Ruderman 1997). Nonlinear MHD calcula- 
tions and further investigations of neutron star physics will be 
needed to resolve these issues fully. 
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